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ABSTRACT 

The magnetic fields and energy flows in an astronomical jet described by our earlier 
model are calculated in detail. Though the field distribution varies with the external 
pressure function p(z), it depends only weakly on the other boundary conditions. 
Individual fieldlines were plotted; the lines become nearly vertical at the bottom and 
are twisted at the top. An animation of a fieldline's motion was made, which shows 
the line being wound up by the accretion disc's differential rotation and rising as a 
result of this. The distribution of Poynting flux within the jet indicates that much 
of the energy flows up the jet from the inside of the accretion disc but a substantial 
fraction flows back down to the outside. 



1 INTRODUCTION 



Collimated radiating jets are a feature of many astronomi- 
cal objects. They have been observed in radio galaxies (see 
Hargrave and Ryle (1974), Carilli and Barthel (1996), Gizani 
and Leahy (2003), Begelman, Blandford and Rees (1984)), 
quasars and young stellar objects. They may also be present 
in gamma ray bursts (see Uzdensky and MacFadyen (2006)), 
dying stars and microquasars (see Mirabel and Rodriguez 
(1999)). In all these objects, the jet emerges along the axis 
of an accretion disc. As the form of the highly relativistic jets 
in active galactic nuclei and the much slower ones around 
stars is very similar, it is likely that relativistic effects are 
not essential for the formation of jets. The presence of strong 
magnetic fields in radio galaxies was deduced from their syn- 
chrotron radiation by Baade (1956), but the similarity of the 
observed structures suggests that magnetism plays a large 
role in the formation of all jets. In particular, recent obser- 
vations of maser emission confirm that magnetic fields are 
important in collimating stellar jets, i.e. Vlemmings et al. 
(2006). 

Following earlier numerical work on the Grad-Shafranov 
equation by Barnes and Sturrock (1972) and Sturrock and 
Barnes (1972), Lovelace (1976) obtained jet-like magnetic 
structures by considering the winding up of a pre-existing 
uniform field. However, collimation by a pre-existing field 
side-steps the problem of why the field is collimated to start 
with. In Lynden-Bell (1996) (paper II), the growing mag- 
netic towers picture was developed in which a magnetic 
field twisted by an accretion disc was confined by an am- 
bient coronal pressure. This provided a crude picture as to 
why the field was collimated. In later work (Lynden-Bell 
(2006) - paper IV) a non-relativistic model was developed 
which involved a force-free magnetic field being twisted up 
by a swirling accretion disc to give a jet-like collimated mag- 
netic tower, now in a stratified atmosphere. Expressions for 



the magnetic fields were given; however, the fields were only 
illustrated in the idealised "dunce's cap" configuration for 
which the solutions were explicit and exact. In this paper the 
more general implicit expressions will be numerically eval- 
uated to depict in detail the fields and energy flows within 
the jet predicted by the model under various boundary con- 
ditions. 

Kato et al. (2004) have computed the starting of 
such jets in a magnetic accretion disc. Li et al. (2006) 
extended these ideas to long jets that are not force-free 
initially but become so as they evolve. Nakamura, Li and 
Li (2006) performed numerical calculations of jets within 
an isothermally stratified atmosphere up to heights that 
are twice the diameter of the disc; our analytical solutions 
extend to far greater heights. Discussions of jet stability 
are given by Nakamura and Meier (2004) and Nakamura, 
Li and Li (2007). Lovelace et al. (2002) computed solutions 
in a toroidal box and found a strongly collimated Poynting 
flux as predicted by the growing towers picture. However, 
their calculations were limited by their box. Romanova 
et al. (1998) and Ustyugova et al. (1999) computed the 
evolution of magnetic fields above discs, with the fields 
not force-free. Though Ustyugova et al. considered fields 
initially of split monopole type and Romanova et al. treated 
a more complicated field structure, both obtained magnet- 
ically dominated collimated outflows with the matter flow 
everywhere away from the disc. As discussed later, we find 
that there is also a smaller Poynting flux back down onto 
the outer feet of the fieldlines, which is due to the fieldlines' 
tension pulling on the disc. There are many simulations of 
highly relativistic black hole environments, such as Gammie 
et al. (2003), McKinney (2005) and McKinney and Narayan 
(2007a/b). 
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2 THEORETICAL BASIS 

A brief summary of the theory in paper IV will be provided, 
as the investigation of the fields which follows is based on the 
model developed in this paper. Cylindrical polar coordinates 
are used. The magnetohydrodynamic equations are assumed 
in mixed Gaussian units. 



2.1 Introducing the Problem to be Solved 

The field is located within a magnetically dominated cavity 
in the external medium. The displacement current is ne- 
glected, so that V x B = 47rj. Though the inside of the 
cavity has negligible gas pressure, it is assumed to be a per- 
fect electrical conductor, so that the medium is frozen onto 
the magnetic fieldlines. Due to the negligible gas pressure 
the magnetic force has nothing to oppose it so the field in 
the cavity is force free (j x B = 0). Hence a solution is de- 
fined by the following boundary conditions: the normal field 
B n at the base, the twist on each fieldline <£, and the gas 
pressure on the outside of the cavity p(z). As the equations 
for the magnetic field do not involve time derivatives, the 
time dependent solution can be found as a series of static 
solutions with differerent twists of the fieldlines at different 
times t. 

Within a radius R at z = 0, a total magnetic flux 
P(P, 0) rises out of the accretion disc. The fieldlines cor- 
responding to this P turn over and return to the disc at an 
outer radius; axial symmetry means they define a surface of 
constant P(P, z) as shown in figure (|4]), where P(P, z) gives 
the flux through a circle of radius R at height z. The disc's 
rotation with angular velocity £} d (R) produces a total twist 
on the fieldlines of $(P) = [Q d (Ri(P)) - Q d (R (P))]t = 
Q(P)t. Here Ri(P) and R (P) are the radii of the inner 
and outer footpoints of the field line on the flux surface 
P(R,z) = P ^constant. At the boundary of the magnetic 
cavity, an external gas pressure p(z) is specified which must 
be balanced by the internal magnetic pressure. 

As V • B = the magnetic field may be written in terms 
of the flux function P(R,z) as (see Mestel (1999)): 

B=ivPxV^ + 5^ (1) 

where <j> is the azimuthal coordinate. With the force free 
condition 



jxB = 0^VxB = 47rj = f3 f B 



(2) 



this implies that taking the poloidal component of the curl 
of (pQ) results in 
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where /3(P) is constant along each fieldline. The toroidal 
component leads to the differential equation for the flux 
function P: 
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]_dP\ 



+ 



OR \RdRj dz 2 



-P(P)(3'{P) 



(4) 



where the magnetic pressure must be B 2 /8tt = p(z) on the 
surface of the magnetic cavity defined by R = Rm(z). This 
equation is nonlinear and is to be solved within an unknown 
surface with an unknown function /3(P) (which must result 
in the specified values for <E>(P)). 



2.2 Variational Method and its Results for the 
Cavity Shape 

The easiest way to overcome these difficulties is to minimize 
the energy W given by 



W 



-ass 



Bj + B 2 R + B 2 
8tt 



RdRd(j) + p(z)A dz (5) 



using a trial function, apply the boundary conditions and 
solve the resulting variational equations as is done in paper 
IV. As the radial flux is constant while the magnetic cavity 
grows to a great height, Br is neglected. Using a trial func- 
tion for the toroidal flux such that the twist on fieldlines 
is the geometric mean between the case with all the twist 
at the top and that with the twist uniformly distributed in 
height, the variational equation was solved generally and the 
following results were obtained: 



ttRI = IP m / x /27rp(z) 



Z) = J 2 P m <S> 2 (P m )/(167Tl) 



(6) 

(7) 



where $ 2 = ($(P m )/P m ) $(P)dP, P m (z) is the maxi- 
mum flux P at a height z and F is the maximum P at z — 0. 
J and I are defined as J = (P|) /SJ and I = (B 2 ) /\B z \ 2 
where () indicates averaging over a horizontal plane, B^ = 
R- 1 f* m B^dR and \B Z \ = A~ x f* m \B z \2nRdR; however, 
we approximate J = I = 1 for our calculations (for exact 
values see paper IV). 

The shape of the magnetic cavity can be found by solv- 
ing equations © and for the cavity radius R m (z). How- 
ever, this only results in valid solutions when the pressure 
falls off more slowly than for n — 4, where p(z) = a(z + b)~ n . 
(paper IV also develops a ram pressure confined model which 
is valid for pressure falling faster than n — 4.) Plots of 
the resulting cavity shapes show that as the jet rises it be- 
comes more collimated; constant pressure gives conical cav- 
ities whose cross section decreases with height, deceasing 
pressure gives centrally bulged cavities. 

The two models we will use for the relative angular ve- 
locity of the footpoints H(P) are 

i) the simple model H(P) = Ho(l — P/F), where Qo is the 
maximum of £1 

ii) the more physically realistic dipole model H(P) = Ho(l — 
(P/P) 3//2 ) - which corresponds to a uniformly magnetised 
rotating body surrounded by an accretion disc. It arises from 
a flux distribution P(P,0) = FR 2 /l 2 for R < I (the "inside 
of the disc") and P(P,0) = Fl/R for R > I (the "outside 
of the disc"), combined with a footpoint angular velocity 
Q d = Q for R < I and Q d = Q (P//)" 3/2 for R>1 



2.3 Magnetic Field in the Cavity 

The field within the cavity is found by writing 

P(R,z)=P m (z)f(X) (8) 

for each height, where A = R 2 /R 2 n (z). The resulting con- 
ditions on /(A) (/(0) = /(l) = 0, /max = 1) are very re- 
strictive, so that in the tall tower approximation one can 
neglect variation of / with height. Additionally neglecting 
second derivatives and squares of first derivatives of R m (z) 
and P m (z), equation Q becomes: 
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4P m d 2 /_ gg , 



(9) 



which gives /3 'ft' as a product of a function of z and A. 
However, it is also a function of P, which is such a product. 
These facts imply that it is a power law. Now / Br is zero 
on axis so /3(P) is zero when P is zero. It follows that since 
/3/3' is a power law, /3 is as well, so we can write (3 — C\P V . 
Inserting this in equation J4} and separating variables gives: 



Rrn — C'2 Prr 



dA 2 



(10) 
(11) 



with C 2 = ClCll A. 

The second equation can be used to find an expression 
for v. Multiplying it through by —//A, integrating by parts, 
and multiplying by P^/A 2 it can be shown that this equa- 
tion becomes 



(*2> = "<*S> 



(12) 



By expanding a slice through the jet as in Lynden-Bell 
(2003) with Br = 0, a different expression for the ratio 
\B 2 Z ) I (Bl) can be derived: 



V(2 -a) 



where at each height 



s{z) = 



1 



top 



A(z)(dp/dz)dz 



A(z)p(z) 

Thus v becomes a weak function of z 
u{z) = l/{2-s{z)) 



(13) 



(14) 



(15) 



so that the weak variation of the profile with height is de- 
termined, as a profile /(A) for every height can be obtained 
from CD) . 

When the external pressure is constant, / can be found 
exactly (as v = 1/2). Integrating equation ([TT]) leads to the 
solution 



/(A) = -eAlogA 



(16) 



When the pressure falls off with height, it is useful to 
obtain approximate analytic solutions for /(A) to increase 
calculation speed. 

In the cases plotted here, the value of v{z) is mostly < 1. 
A good approximation for v(z) < 1, which we will use for 
falling pressure throughout this paper, is /(A) = g(X)/g(Xi) 
where 



3(A) = 



A(l-A a ) 
a(l + a 2 \ a ) 



(17) 



with CL2 determined by the condition that / has a maximum 
at A = Ai 

a 2 = (A?(l + a)-l)/(A?((l-a)-A?)) 

and witha(z) = 2v(z) - 1 and Ai(z) = (1 + 0.07073a(z))/e. 
Comparing this solution to exact numerical computations 
shows that it is gives good results for v < 3. 



3 ELECTROMAGNETIC FIELDS 
3.1 General Method 

The jet fields must be found numerically in all but the sim- 
plest cases, as the fields are described by both implicit and 
differential equations. The calculations were performed us- 
ing the Mathematica program. The errors, unless otherwise 
indicated, are negligible compared with the magnitude of 
the results. 

The different boundary conditions which define the 
model are specified by the parameters £, Oo, F as well as, as- 
suming a power pressure law, n, a, b where p(z) = a(z+b)~ n . 
As similar jets occur in a variety of objects with very dif- 
ferent sizes and properties, we decided not to attempt to 
choose physically realistic values for the parameters; instead 
we use the following default values to investigate the fields: 
We choose units of time t such that 2ty corresponds to one 
turn at angular velocity Qo and measure magnetic flux in 
units of F. The lengthscale was then chosen so that a = 1. 
With b = 0.1 this corresponds to a radius of the disc of or- 
der b for the n — 3 case. Most plots are either for n — 3 
(which is a realistic value) or for constant pressure. Values 
of t = 30 and b = 0.1 were used; these give a collimated 
tall tower where the approximations are valid. We use the 
simple model H(P) unless mentioned otherwise. 



3.2 Flux Surfaces 

As the fieldlines rise up from and return to the disc, they 
define a tube-like surface of constant flux P. The equation 
giving P is 



P{R, 



Pm(z)f(X) = R 2 m (z)^27rp(z)f(R 2 /R 2 m (z))/I(18) 



The cross sections of the surfaces of constant P on which the 
lines of force lie are plotted in figure (pQ) for both constant 
pressure and for a pressure which falls as n — 3. The bound- 
ary of the magnetic cavity is also indicated in the diagrams. 



3.2.1 Results 

For the n — 3 and n — external pressure distributions 
the constant flux contours / poloidal fields conform to the 
shape of the corresponding magnetic cavity. At the centre of 
the jet, the contour lines are nearly vertical; they are closely 
spaced at the base. 



3.2.2 Discussion 

In the constant pressure case, it is remarkable how simi- 
lar the form of the contours and the cavity are to those 
in the "dunce's cap" model evaluated analytically in paper 
IV. While the pressure in the "dunce's cap" model is also 
constant, H(P) on the disc is different to the simple model 
used here. This result suggests that the field structure and 
the shape of the jet are determined mainly by the external 
pressure as the disc's rotation profile only has a secondary 
influence on the jet's structure. 
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Figure 1. Cross sections of surfaces of constant flux, which also 
give the lines of force with the azimuthal component suppressed. 
On left, for p oc [z + 0.1) -3 at t = 30; on right, for constant p at 
t = 90. Both plots for are for the simple model relative rotation of 
footpoints. (Note: in all plots lighter shades correspond to larger 
values.) 



3.3 Magnetic Field Components 

Substituting equation (|8} into the expressions for the fields 
given by equation Q leads to 

B z (R,z) = ^f'(X) (19) 

B ^z) = ^( P rn(z)f(\)y^ (20) 

with P m (z) = 7vRl l (z)^27rp(z)/I and A = R 2 /R 2 m {z). 

\B Z \ and \B<f,\ were calculated for the n = 3 pressure 
law and are shown in figure J2}. B 2 — P| + B 2 and the 
twist | | / 1 | were plotted for both constant pressure and 
the n = 3 pressure law in figures (|2} and J3} respectively 
(n = plots were made at t = 50 so that both cavities have 
nearly equal height). The plots are linear; to preserve detail 
we cut them off at large field strengths (n = 3 cutoffs 270, 
10, 15 for P 2 , l-B^I, \B Z \). To investigate the dependence of 
the fields on Q(P), B 2 was also found for the dipole model 
Q(P). It is depicted in figure J2}. 

It is also useful to determine how the fields change when 
the values for t, Qo, F, a and b are varied. For the n — 3 case 



the parameters were individually changed from their default 
values and the results noted. 

3.3.1 Results 

When the pressure falls as n — 3 it can be seen that B z is 
largest at the base and core of the jet. At a radius Ri(z) 
which corresponds to A i at each height, B z goes to zero. B$ 
is large at R\ and at the base; there is almost no azimuthal 
field at the centre and at the edges of the cavity. The twist is 
highest at R\ , and is negligible on axis and near the bound- 
ary. B 2 is maximal at the bottom of the jet and on the upper 
axis. 

The constant pressure jet has the highest B 2 field 
strength on axis. The field's twist is zero near the axis and 
the boundary. 

The plot of B 2 for the dipole rotation model is very sim- 
ilar to the simple model plot; the most noticeable difference 
is that the dipole jet is taller. 

The height and collimation of the magnetic cavities in- 
creases with time t, though the distribution of fields does 
not change much. Varying the other parameters Qo, F and 
a changes the height of the cavity h c & v at a given time (in- 
creasing Qo or F increases the height, increasing a decreases 
it); however, this only affects the field distribution substan- 
tially if the cavity is thus in a different pressure regime. So 
for p(z) = a(z + 6) -3 , if h cav becomes < 6, the fields become 
as in the n — case, if h cav « 6, the fields are an intermedi- 
ate case, and when h cav > b we retain the fields as in figure 

©• 

3.3.2 Discussion 

One clearly expects strong magnetic fields at the base be- 
cause all the fieldlines rise up from and return to the disc. 
When the outside pressure is ~ constant, theory predicts 
pinching of the central axial fields by the toroidal fields, 
which results in a high field strength on axis, as is seen in 
figure ©. In the n — 3 jet this pinching is only seen at the 
top, where s(z) becomes small and v(z) approaches 1/2 as 
in the constant pressure case. The strong central fields may 
cause the emission of most of the radiation in some of the 
observed jets. 

The behaviour of B z , B^ and the twist at Ri is due to 
the fieldlines turning over at this radius with a large propor- 
tion of their total twist at their tops. 

The observation that the form of the fields does not 
depend very strongly on H(P) or the other boundary condi- 
tions lets us deduce that the plots made are quite generally 
valid. Ram pressure confined jet fields (see paper IV) are 
probably like the n — 3 fields, as the cavity shapes are simi- 
lar. The increased cavity height for the dipole model at the 
same t is due to the larger differential rotation of the dipole 
model for all values of P. 

3.4 Individual Fieldlines 

In order to see the shapes of individual fieldlines, 3D plots 
were necessary. The field lines are given by 6.R/Br = 
R&cjy/Bcf) — dz/B z , or using the expressions for the com- 
ponents of B: 
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Figure 2. First four from left: plots of \B Z \, \B (f) \, \B (f) \/\B\ and B 2 for simple model Q(P). For comparison, rightmost plot shows B 2 
for dipole model disc rotation. All plots are at t = 30 with p oc (z + 0.1) -3 



dR 



dz 



dP/dz 0(P) dP/dR 

Thus a differential equation for the fieldline is: 

dz _ n 9P/dR _ , jir_ 



(21) 



(22) 



The rightmost expression is a function only of z for a given 
fieldline, i.e. a given constant value of P. The expression is 
double valued: the fieldline rises from and then returns to the 
disc, so that X(z) is double valued. Two different functions 
for X(z) were obtained by solving f(X,z) = P/P m (z) for 
X(z) subject to either Xi n (z) < Ai or X ou t(z) > Ai. The 
function z{<j>) for all points on the fieldline was found by 
solving equation ([22]) numerically using X- lxl (z) up until the 
point where z reached a value corresponding to the top of 
the fieldline, then integrating back down to the base using 
Aout(^). The corresponding values of R((/)) were found using 
R 2 = R 2 rn (z(0))X(z(0)). 

3D plots of the P = 0.3F fieldline (for n = 3 and simple 
model n(P)) were thus constructed at three different times. 
They are shown in in figure JD along with the constant P 
surface and an outline of the magnetic cavity. 

In order to be able to visualize the time development 
of the fieldlines, a 3D animation of the motion of the P = 
0.3P fieldline was made for 30 < t < 80. The outer foot 
point of the fieldline was set to rotate with the small angular 



velocity = 0.07. The animation is provided on the website 
"www. ast . cam .ac.uk/ research / repository /bfieldline . ht ml" . 



34.1 Results 

The fieldlines are twisted most at the top; the outer part of 
the fieldline has less twist on it than the inner part. 

The video depicts the fieldline being wound up and ris- 
ing rapidly. At the bottom of the cavity the fieldline's twist 
decreases, so that eventually the line is nearly vertical. 



3.4-2 Discussion 

The increasing twist near the top is due to the trial function 
used. In the video it can be seen that the fieldline is dragged 
around by the disc at the inside foot point and held back by 
the outer foot point. The fieldline's tension thus applies a 
torque and transfers energy from the inner to the outer part 
of the disc. 

The straightening of the fieldlines near the base means 
that the fields become nearly antiparallel, which could lead 
to magnetic reconnection if the perfect conductivity breaks 
down. 
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3.5 Poynting Flux 

As there is no E along B due to the assumption of perfect 
conductivity, and as the velocity of a fieldline u is the cE x 
~B/B 2 drift velocity, the electric field is given by E = — u x 
B/c and u • B = 0. The azimuthal angle of an intersection 
of the fieldlines with a z ^constant plane is shown in paper 
IV to be 

= * + n d (Ro)t + n(P)tq(z, P) (23) 

where Qd(Ro) is the angular velocity of the outer foot point 
of the field line, \£ is the initial azimuth at which the 
fieldline intersected the disc at the outer foot point, and 
q= J ^(dlog A/dlog/)dz/ J line (d log A/d log f)dz is the frac- 
tion of the total twist on the field line P which occurs by 
height z (integrating up from the outer foot point). In the 
same paper, the velocity of the intersection is derived in 
terms of the components of the fieldline's velocity u: 

R<j> = U(f> (l + Bl/B 2 z ) + u R (B R B^/B 2 z ) (24) 

Applying Faraday's law to the flux through a circle about 
the axis in the z ^constant plane gives 

= -P/c = -(u z B R - u R B z )/c (25) 

We solve these two equations and u • B = for u; neglecting 
Br we obtain 

—B^BzRcj) -u z B z -P 

Uz= \B\* 5 U * = -BT ; 
E follows directly from cE = B x u, giving 

cE R = -R(j)B z ; cE (f) = -P; cE z = -PB^/ B z (26) 

It is interesting to investigate the distribution of the Poynt- 
ing flux component S z . As S = ^E x B 

S oc \B\ 2 u 

so that 

S z oc —B(t>B z R(j) (27) 

q was approximated as Aq/At by calculating Aq at times 
separated by At = 0.03 (using either Ai n or A ou t to overcome 
the double- valued-ness of q) . Values of t = 50 and b = 4 were 
used to calculate S. H(P) was given by the dipole model as it 
has a well defined Qd • It was verified that the changes did not 
affect the field distributions substantially. As calculations of 
Rcj) were very time consuming, B z R(j) was evaluated on a 
70 x 70 evenly spaced grid across the width and height of 
the jet; cubic interpolation was then applied to obtain values 
at other points. The vertical Poynting flux was thus plotted 
in figure JS}. 

3.5.1 Results 

The z component of the Poynting flux is large and positive 
just above the inside of the accretion disc. Above the outer 
part of the disc there is a large downward flow. 

3.5.2 Discussion 

Energy flows up from the inside of the disc because the disc 
drags the fieldlines around and thus does work on the field. 




Figure 5. Plot of RS Z for the dipole model at time t = 50 with 

p oc (z + 4) -3 

The flow spreads to the top and sides of the jet, where it 
enables the expansion against the external pressure. Some of 
the energy returns to the outer part of the disc (as the fields 
pull on the disc there) ; the fraction this constitutes was cal- 
culated as f*™ (0) RS z dR/ f* 1 RS z dR = 0.47. It thus ap- 
pears that jets are quite an efficient mechanism for energy 
transfer from the inside to the outside of the accretion disc. 
The energy flow into the jet can also be calculated by consid- 
ering torques on the disc. The couple GdP — — 2nRdP 
carried by the tube of flux between P and P + dP is con- 
stant along the tube. This couple inevitably does work on 
the rotating outer feet of the field lines, giving a power flow 
G{P)dP Q(R (P)) down from the jet into the outside of 
the disc. This differs from the results of many simulations 
(Ustyugova et al. (1999), Gammie et al. (2003), McKinney 
(2005), McKinney and Narayan (2007a/b) and Narayan et 
al. (2007) on winds) because they have no magnetic flux re- 
turning to the disc and so do not have the downward Poynt- 
ing flux which results from this. Romanova et al. (1998) do 
not have a force- free magnet osphere. The disc does work 
on the field lines at the inner feet, giving a power input of 
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4 CONCLUSIONS 

1. The B field is strongest at the base of the magnetic cavity 
and on axis. The field structure is only weakly affected by the 
detailed profile of the shear in the accretion disc's rotation. 

2. The fieldlines are twisted more at their tops and at smaller 
radii. At the cavity base the fieldlines become straight and 
vertical. 

3. In the case investigated the energy rises from the inside of 
the disc; somewhat less than half returns to the outer part 
of the disc, the remainder flows up the jet. 

The published numerical calculations of others referred 
to above do not extend to great heights due to the difficulties 
with numerically induced resistivity. Although in our work 
only the "dunce's cap" solution is truly exact, the solutions 
in this paper are still basically analytic and do not suffer 
this inconvenience. In the few cases such as Nakamura, Li 
and Li (2006) where similar problems are solved, the results 
agree with the start of our solutions. 
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G(P)dP Q(Ri(P)). The difference between the power input 
and output goes into expanding the magnetic cavity against 
the pressure and increasing the magnetic energy. Integrating 
over all flux tubes, the net power put into the expanding jet 
is/ F G(P)O(P)dP. 

As —P/c = E<j> is small compared to the other E com- 
ponents but EJ_B, E points close to normally out of the P 
surface at each point. 



